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Abstract 

It is clear that the solar corona is being heated and that coronal magnetic fields undergo 
reconnection all the time. Here we attempt to show that these two facts are in fact related i.e. 
coronal reconnection generates heat. This attempt must address the fact that topological change 
of field lines does not automatically generate heat. We present one case of flux emergence where 
we have measured the rate of coronal magnetic reconnection and the rate of energy dissipation 
in the corona. The ratio of these two, P/$, is a current comparable to the amount of current 
expected to flow along the boundary separating the emerged flux from the pre-existing flux 
overlying it. We can generalize this relation to the overall corona in quiet Sun or in active regions. 
Doing so yields estimates for the contribution to corona heating from magnetic reconnection. 
These estimated rates are comparable to the amount required to maintain the corona at its 
observed temperature. 


1 Introduction 

There is still little consensus on what mechanism can be credited with supplying heat to the Sun’s 
corona. Among the most frequently invoked candidates are dissipation of waves and magnetic 
reconnection. Both processes are known to occur, but their relative contributions to heating has 
yet to be definitively quantified observationally. 

There is extremely strong evidence that magnetic reconnection is occurring throughout the 
corona at some rate. The coronal field is connected to photospheric flux concentrations which are, 
in all the best observations, surrounded by photosphere unconnected to the coronal field, if not 
entirely unmagnetized. These flux concentrations move about, apparently at random, under the 
influence of granular and super-granular flows. If the coronal field lines remained anchored to the 
same pair of footpoints over days or weeks, the coronal magnetic field would appear extremely 
tangled and complex. The coronal field outlined in EUV images shows little sign of such tangling 
— in fact, it appears smooth enough to have been “combed”. While it is still possible that complex 
tangling occurs at length scales below our present resolution [1], footpoint motions occur over all 
length scales and presumably so should the tangling. There is little evidence for it on the the largest 
scales, which appear increasingly well fit by potential fields over time [2]. This fact gives a clear 
indication that coronal field lines are constantly being reconnected: uprooted from one footpoint 
and reattached to another. 

Based on this reasoning we propose that magnetic field lines in a given portion of the corona 
are undergoing topological change at some rate T. Two questions raised by this proposition are, 
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firstly, how is topological change related to heating?, and secondly, what fraction of overall coronal 
heating can be attributed to topological change? We address these two questions in the following. 

2 Heating from magnetic reconnection 

To illustrate how the topological change of magnetic field lines could result in plasma heating we 
consider a simplified, two-dimensional model of flux emergence sketched in Fig. [1] In the initial 
state, Fig. [ljr, a small bipole, P2-N2, sits beneath a larger, existing bipole Pl-Nl, and the coronal 
magnetic field is potential: V x B = 0. The next three panels, [Do |T]d, show a state in which the 
small bipole has emerged further (P2 and N2 have separated and their fluxes have increased) while 
the outer bipole has not changed. 






Figure 1: Stages, progressing clockwise from upper left, in reconnection following the emergence of 
the bipole beneath existing bipole. Field lines are blue curves originating in sources PI, P2, Nl, 
and N2, located in the photospheric boundary (black horizontal line), (a) Shows the pre-emergence 
state, (b)-(d) show stages after emergence of bipole P2-N2 — (b)-(d) have the same photospheric 
field, (b) is the state before any reconnection, (c) is after some reconnection and (d) is the state 
after complete reconnection. Red lines show the separatrices, magenta curves are the current sheet. 
The shaded region shows the flux which is reconnected to produce state (d) from state (b). Dashed 
curves show a surface though which flux connecting P2 to Nl can be computed. 

In Fig. [Hi, where the coronal field is once again potential, the flux interconnecting old and new 
poles, P2 and Nl, has increased. The interconnecting flux is computed through a surface, depicted 
by a dashed curve, which extends from the coronal X-point to a point A located on the boundary 
somewhere between P2 and Nl. Between Figs. Q}t andQJl this flux has increased by the addition of 
flux shown shaded in [Dl - According to Faraday’s law this additional flux must have arisen through 
an electric fiek@ 

6 = E-dl , (1) 

1 We use cgs Gaussian units with the exception of electric field and current. In the interest of brevity we use for 
those, E = cE cga and / = I C gs/c, whose units are G • cm/s and G • cm respectively. 
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along the perimeter C of the surface. If we assume no interconnecting field lines moved upward 
through the boundary at point A (such motion occurred only within the region between P2 and 
N2, as flux emergence) then there is no electric field along the horizontal leg at A. Instead, the 
entire electric field in integral (JT]) occurs at the coronal X-point — this is the reconnection electric 
field. 

If the coronal field had remained entirely potential throughout the emergence the reconnection 
would have produced no coronal heating at all — there would have been topological change in the 
absence of heating. The rate of electromagnetic work on the plasma is 

P = I<f> , (2) 

where I is the current flowing along perimeter C. Without coronal current there is no electromag¬ 
netic work, and therefore no heating from the topological change. 

Heating occurs only when the magnetic reconnection is slow enough to permit the accumulation 
of current at the X-point. FigureQ}:* shows the case where no reconnection at all has occurred during 
the flux emergence. The flux interconnecting P2 and N1 is therefore the same as before emergence 
and there is a current sheet separating the newly emerged from the pre-existing flux [3j. (The state 
with a current sheet is the one with the lowest magnetic energy subject to the single constraint on 
the interconnecting flux [4].) Any electric field within the current sheet will increase the P2-N1 
flux, taking Fig. [TJo to Fig. [TJc. In this process the magnetic energy will decrease by doing work on 
the plasma, endowing it with either heat or bulk kinetic energy. The rate of electromagnetic work 
is given by Eq. ([2]) where I is the current in the coronal current sheet at which the reconnection 
occurs. 

3 A case study: flux emergence in AR 11112 

The foregoing simplified example illustrates how flux emergence makes topological change clearly 
identifiable, and even quantifiable. Longcope et al. [5] measured such a change for one active 
region (AR) emerging in the vicinity of an existing AR. This allowed them to compute the rate of 
topological change $. Unfortunately, the limited EUV and X-ray data available made it difficult 
to compute a heating power P for this case (TRACE was observing at high cadence, but only 
in 171 A, and Yohkoh SXT had multi-filter partial frame images only during two intervals). The 
observation was thus not ideally suited to understanding the relation between reconnection and 
heating captured in Eq. Q. 

Recently Tarr et al. [6] used the AIA and HMI instruments on SDO [i7j 11] to observe a magnetic 
bipole emerging within an existing AR. Figure [2] shows 21 lA imaged from AIA of the emerging 
bipole (top row) along with radial field maps from HMI (bottom row). The EUV images clearly 
show a dome of flux anchored in the newly-emerged positive polarity, which we call P2 by analogy 
with Fig. [I] The newly emerged negative polarity has moved to the southeast, but the dome clearly 
includes negative polarity to the west: old flux, outlined in blue in the bottom row of Fig. [2], and 
hereafter called Nl. The distinction between old flux (Nl) and new flux (N2) is made by tracking 
the evolution of the magnetograms, first automatically and then adjusting manually [9, 6]. This 
introduces the largest source of error into our calculation. One indication of its magnitude is that 

2 The 21lA filter of AIA is primarily sensitive to Fe xiv, whose peak formation temperature is T ~ 2 x 10 6 K, in 
ionization equilibrium [7]. 
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the total signed flux attributed to newly emerged flux remains constant within 10% of T 2 over the 
primary emergence period — Oct. 15. 
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Figure 2: The emergence of a bipole into the negative polarity of AR 11112. The top row shows 
AIA 211A images from three times during the emergence, plotted with a logarithmic grey-scale. 
The bottom row shows HMI maps of B z (x, y ) from the same times (on a grey-scale capped at ±500 
G). The green curves outline the dome from the 211 A images, and the blue curves outline the old 
negative polarity, Nl, to which reconnection occurs. 

The perimeter of the 211 A dome was traced (manually) at one hour intervals and mapped onto 
the HMI radial field maps (green curves). The total unsigned flux within this boundary remains 
constant to with 10% of ^ 2 ; suggesting that the boundary accurately identifies a closed magnetic 
system [6j. The region of old negative flux (Nl) lying within the dome perimeter was taken to be 
reconnected flux, analogous to the P2-N1 flux in Fig. [I] The integral of B z over the overlapping 
region gives the reconnected & r (t) shown as a red curve in in Fig. [3l A crude linear fit to this 
(dashed line) rises at <h r ~ 2.9 x 10 1 ' 5 Mx/s (Tarr et al. J6| fit with a line starting later and thus 
found a slightly larger slope). 

This value can be compared to that from a potential held extrapolation from the same radial 
held. A magnetic charge topology (MCT) extrapolation replaces each region by a point source in 
order to tally the total flux connecting every pair of regions [flO] , The hux connecting the emerged 
positive (P2) to the existing negative (Nl), outlined by blue contours in Fig. [2], called \t), is 
shown as a green curve in Fig. [3j and rises at ~ 4.7 x 10 15 Mx/s (dashed line). This can be 
taken as an upper bound on the reconnected hux, in the same way that the potential held in Fig. 
QJ1 included more inter-connecting hux than the states with current, Fig. [T}o or[l]x 

The emerging region was also observed by the X-ray Telescope (XRT) on Hinode [TTj. A ratio 
of images from Ti-poly and Al-mesh hlters is used to derive the temperature and emission measure 
of the plasma within the dome PEI El- This is then used with the total radiative loss function 
from Klimchuk and Cargill m to derive the power, P r . radiated by the dome plasma, shown as 
a black curve in Fig. [4l To generate this power entirely by the reconnection at the observed rate 
<h r ~ 2.9 x 10 15 , would require current / = P r / ( h r , in accordance with Eq. Q. This value, found 
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Figure 3: The fluxes from the emerging bipole over time. The blue curve shows the total flux in 
the positive polarity P2 in units of 10 21 Mx. Diamonds mark the times of the panels in Fig. [2j The 
red curve is T r (i) found from integrating the portion of the old negative polarity within the dome. 
The green curve is the amount of reconnected flux in a potential field, Dashed lines are 

linear fits, both intersecting zero at 11 : 15 on Oct. 14. 


by reading the black curve from the right axis, rises to over I~3x 10 11 Amps over the 50 hours 
of emergence. 

In a simplistic view of this process, the emerging flux reconnects to its surroundings in order to 
reduce its magnetic energy, and thereby approach a potential field. The reconnection is therefore 
driven by the difference in fluxes between the actual interconnecting flux and that in a potential 
field 

AT = T< v) - T r , (3) 

which is the separation between between the green and red curves in Fig. [3l This flux difference is 
plotted in violet in Fig. [2 In order to plot it with P r . the flux difference AT is multiplied by the 
scaling factor £ = 1.6 x 10 5 G cm/s, chosen to bring the two curves into approximate alignment in 
the plot. (A possible interpretation of this empirical scaling factor is suggested below.) 

If the reconnection were occurring across one or more current sheets created only in response to 
the flux discrepancy AT then we would expect a relation between AT and /. The simplest relation 
would be a linear one 

_ AT 
““ ~C 

where C is formally a self-inductance, although the analysis so far has no current path with which 
to associate it (one possibility is suggested below). Equating this current with the value, P r / T r , 
derived from the power, yields an expression for the inductance 


C 


T r AT T r 

_ _ 

Pr ~ C 


1.8 x 10 10 cm , 


(5) 


after using the values quoted above. 
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Figure 4: The radiated power P r from the dome plasma computed from Hinode XRT filter ratios. 
The power is a black curve read in units of 10 25 erg/s against the left axis. The same curve is scaled 
by 4> r ~ 2.9 x 10 15 to yield a current, read from the inner right axis in units of 10 11 Amps. The 
flux difference, defined by Eq. ([3]), is plotted in violet after being scaled by the factor £ = 1.6 x 10 5 
G cm/s. It can be read from the outer right axis in units of 10 21 Mx. 


A linear relationship like Eq. (J4]) obtains in the simple two-dimensional model of Fig. |T[ the 
current in the single, global sheet is proportional to the amount of grey, unreconnected flux under¬ 
neath it HU. In a three-dimensional model with isolated sources, the topological role of the X-point 
is played by a separator, and a sheet will form there carrying a current proportional the the flux 
difference A<h [ 15] . In that model a separator loop of length L has self-inductance is C = 47tL 
times a logarithmic factor typically of order unity. 

The MCT model of the emerging bipole, shown in Fig. [51 includes a positive coronal null point 
(triangle) 15 Mm above the positive polarity P2, whose fan surface forms a separatrix dome. The 
base of this surface is a ring of spines (solid magenta curves) linking the negatives sources to 
negative photospheric nulls, shown as green and cyan As. These nulls are linked to the coronal 
null by separators lying within the separatrix dome. Two of these separators, shown in red in Fig. 
[5] form a loop enclosing the reconnected flux linking P2-N1 (blue field lines). The current sheet 
will form along this loop, in analogy to that in Fig. [TJ The separators are 24 Mm and 27 Mm 
long, creating a loop of total length L = 51 Mm, comparable to the distance inferred from the 
self-inductance in Eq. ©, although it would seem the logarithmic factor in this case is ~ 3.6. 

The rate of observed magnetic reconnection, <f> = 2.9 x 10 15 Mx/sec, is equivalent to an electro¬ 
motive force (EMF) of 29 Megavolts. If this EMF were from a simple loop, then the mean electric 
field, 4>/L ~ 0.6 V/m, would be about one hundred times greater than the Dreicer field. We do not 
believe that such a large electric field could actually be present. Aside from theoretical difficulties 
implied by such an electric field, there is absolutely no evidence of particles having energies even 
close to 29 MeV. This suggest to us that the electric field within the reconnecting current sheet 
is far more complex than a simple closed loop. The EMF might, instead, be built from numerous 
parallel reconnection events transferring reconnected flux through the separator loop. 


6 





Figure 5: The MCT model of the field at 2010 Oct. 15 13:00 (the same field as the middle panel of 
Fig. [2]). The left shows the points source white xs and black +s on top of the HMI radial field map 
(grey scale). Photospheric negative nulls are indicated with cyan A and the coronal positive null 
with a cyan V. Magenta solid curves show the spine lines from the photospheric nulls. These form 
a closed curve in the photosphere which is base of the separatrix dome — the fan surface of the 
coronal null. The solid orange curves show the two spines from the coronal null (one terminates at 
P2 the other connects to a positive source far away). The dashed magenta curves are the bottoms 
of the separatrix surfaces which intersect the dome surface along the the two separators (solid red) 
enclosing the P2-N1 field lines. The right panel shows the same features in a projection, with 
vertical scale exaggerated for clarity. Nulls are green rather than cyan and a selection of P2-N1 
interconnecting field lines are shown in blue. 


This exercise provides one example where the reconnection rate, <f>, can be measured along 
with the heating power, or at least the portion P r lost through radiation. The values found from 
these observations are in conformance with theoretical models in which reconnection produces the 
heating. In particular, the ratio P r /T, corresponds to a current which can be understood as that 
resulting from the emergence and accumulation of magnetic free energy. 

Having related reconnection to heating in a single case, we take up the task, in the next section, 
of applying the relation to a more general setting: heating the entire corona. Topological change 
can be quantified observationally in cases of small flux elements interacting in the quiet Sun. The 
heating from these interactions can be used to predict a mean coronal heating rate from quiet Sun 
reconnection. We apply a similar technique to predict heating from reconnection in a active region. 

4 A statistical view of reconnection heating 

To estimate the net contribution of reconnection to coronal heating we consider a collection of 
photospheric flux concentrations, of mean flux d<b, randomly distributed with areal density n. 
To remain in steady state one concentration must completely reconnect in the time r ex that it 
takes it move into the position of a neighboring concentration. If it moves randomly with typical 
photospheric velocity u p h, this exchange time will be 

i Jm/b 

Tex ~ irr — 

V n V ph v ph 
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where B = n5& is the mean flux density of one sign of concentrations. (At this point we treat 
unipolar and mixed-polarity distributions simultaneously.) The rate of flux transfer into or out of 
domains anchored to a single concentration is therefore 


ou* /— — 

<f> = — ~ v ph V B 5$ . (7) 

7"ex 

According to Eq. (f2j) . flux transfer across a current sheet with typical current / will do electromag¬ 
netic work 

P = ~ u ph /V bI , (8) 

on the plasma occupying the coronal field anchored to that flux element. The heat flux deposited 
into the entire corona is therefore 

F h ~ nP = v ph IB 3/2 d$~ 1/2 , (9) 

per photospheric surface area. When the coronal field evolves quasi-statically the current, I, will 
be independent of velocity so the heat flux scales linearly with velocity. This is natural for a quasi¬ 
static heating process which must release a fixed amount of energy from a given displacement, 
independent of how rapidly the displacement is made. 

4.1 Quiet Sun 

The quiet Sun corona lacks large-scale structure, and therefore contains current sheets only between 
neighboring flux concentrations. These sheets extend a distance comparable to the separation 
between elements, L ~ h” 1//2 . They are not driven by large-scale forces, but rather exist to 
compensate for an unreconnected flux which is some fraction / of the entire element. Complete 
flux exchange is therefore assumed to occur in 1// distinct reconnection events. (The value of / 
will depend on the details of the reconnection process, and is considered a free parameter at this 
time.) The current will therefore scale as 


I ~ ~ fy/B6$ . ( 10 ) 

Substituting this into expression (JHJ) gives the heating power into a single element 

P ~ fv ph BM . (11) 

It is evident that little heating power results when reconnection occurs relatively smoothly, in very 
small steps (i.e. / —>• 0). 

In a survey of quiet Sun magnetograms from MDI Longcope and Parnell pj6] found the isotropized 
Fourier power spectrum and the kurtosis of the magnetograms could be matched by discrete flux el¬ 
ements with mean size <54 = 1.0 x 10 19 Mx distributed with an areal density of n = 7.0 x 10” 19 cm” 2 . 
Taking their random velocity to be u p h = 250 m/s [T7j means they would exchange their flux in 
r ex = 13 hours. Close et al. [IB] used a potential field extrapolation from a sequence of from 
MDI quiet Sun magnetograms to compute the reconnection between flux elements. They found 
a reconnection recycling time only one tenth the value derived above. This faster reconnection 
might be due to a faster motion, or from reconnection sequences more complex than the binary 






case considered here. In any event, we take this as evidence that ours is a conservative estimate 
of the quiet Sun reconnection rate. According to the values from Longcope and Parnell, the mean 
flux density is B = n5& = 7.0 G in typical quiet Sun. 

Longcope et al. [£t9| surveyed X-ray bright points (XBPs) using MDI magnetograms and EIT 
images in two wavelengths. They used the EIT data to derive the emission measure and temperature 
of the XBP from which they computed the total radiated power. The same data can, however, be 
used to compute the power conducted downward from the XBP 


Tcond — 


kT 


(W2) 2 


V 


( 12 ) 


where k = ko T 5//2 is the Spitzer thermal conductivity (ko = 10~ 6 in cgs units). It is probable 
that heat conduction occurs by classical means in events of such small energy flux. The separation 
between the magnetic concentrations is L^ p and the short and long axes of the EIT feature, t\ 
and £ 2 , are used to estimate the coronal volume at V = (47t/ 3)£ 2 ^2- Figure [6] shows the conducted 
power from the 285 XBPs in the survey vs. the mean flux of their photospheric concentrations. 
The dashed line shows Eq. m with / = 0.1, Wph = 250 m/s and B = 7 G, as described above. 
The line lies along the top of the observed distribution, suggesting that our values of velocity, mean 
flux density, or /, are slightly larger than reality: their product may be excessive by about a factor 
of three. 
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Figure 6: Power (conductive) vs. flux for a survey of 285 XPSs from Longcope et al. [19j. Conducted 
power, P cond , from Eq. O, is plotted against the flux from an average of both poles in the bipole. 
The dashed line shows Eq. (ED- 

The average reconnection heat flux to the quiet Sun is found using the current for binary 
interactions, Eq. m, in expression @. The result 

4 QS) ~ fv ph B 2 , (13) 

is, remarkably, independent of element size. The values quoted above yield a heat flux Fh = 
10 5 ergcm -2 s^ 1 , consistent with the heat flux inferred for quiet Sun regions [20]. 
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4.2 Active region heating 

Unlike the quiet Sun, an active region has global currents driven on large scales. A force-free 
magnetic field will have a global current density J = aB. The field is anchored to discrete elements 
but we will assume that current does not flow into any of the photospheric elements. Instead the 
coronal current will be carried along separatrices where it can close along the chromospheric canopy 
between concentrations [21] . The result will be coronal current sheets confined to the peripheries of 
cells. The peripheral sheet associated with a single element will carry all the current which would 

have otherwise entered that cell _ 

<5<I> _ 

I ~ J^- = a5<S> . (14) 

B y ’ 

Flux transfer into or out of the domain anchored to this concentration must occur across this 

current sheet, thereby giving rise to electromagnetic work and heating. Substituting this into Eq. 

(J9j) gives an active region heating rate 

4 AR) ~ aUph^M 1 / 2 . (15) 

The differences between this expression and the quiet Sun heating rate, Eq. (iT3|) . are due to the 
different natures of the currents across which the reconnection occurs: generated globally vs. locally. 
One implicit consequence of this difference is that since the AR current is not driven by a flux 
discrepancy, reconnection will not necessarily reduce it. 

The active region heating is also greater due to the greater density of flux concentrations: 
typically B ~ 100 G in plage. Assuming the same mean concentration size, ciT = 10 19 Mx, with 
the same random velocity, u p h = 250 m/s, and taking a typical twist a rs j 10 10 cnT - 1 mm, 
yields a heat flux R ^ 10 7 ergcm 2 s 1 , comparable to that expected of active regions [241120| . 

5 Discussion 

If heating power P were generated by magnetic field lines undergoing topological changes, i.e. 
reconnection, at a rate 4>, then the ratio of these rates P/<3? would have units of electric current. 
The foregoing showed examples where that relation arose from an actual current in the coronal 
magnetic field. In those examples the reconnection electric field does electromagnetic work on the 
plasma at a rate P = / 4>. This is how magnetic reconnection might heat the solar corona. 

From that basic scenario we derived scaling laws quantifying the heat that could arise as mag¬ 
netic elements move randomly over the photosphere, and coronal reconnection occurs to keep the 
field there from becoming excessively tangled. The process described is, at its root, the same one 
used by Parker to arrive at the reconnection heat flux Fh ~ u p h-B~tan0, for field lines pushed to 
an angle 9 from their potential state [25]. Indeed, the relation found here for quiet Sun heating, 
Eq. (1131) . has the same scaling. In our expression, however, B is the density of small-scale flux 
elements, rather than a vertical magnetic field strength. In point of fact, we tacitly assume the 
concentrations have local field strength B z B in order that they be distinct entities. Notably, 
the local field strength does not enter our estimate of heating, since heating is assumed to occur in 
the corona rather than within the photospheric concentration. 

Expression (I13jh and its interpretation, matches the one derived by Longcope and Kankelborg|26] 
under the consideration of interacting nrixed-polarity magnetic elements. That earlier derivation 
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used an interaction model in which sources disconnected from an overlying background held in 
order to reconnect with one another m- The derivation proposed here assumes only that sources 
exchange flux across current sheets created through their binary interaction. 

We recognize that coronal currents have different origins in ARs than in the quiet Sun: the 
former being globally rather than locally driven. This recognition led us to a reconnection heat 
flux, Eq. (USD, scaling as Fh B 3 ! 2 in ARs. This expression followed, once again, from reconnection 
between small magnetic elements, of a single polarity, composing the AR. Priest et al. [28] treated 
a similar scenario in their flux-tube tectonics model, but found a heat flux scaling similar to Parker’s 
(their equation [53]). They assumed that currents between interacting elements were driven locally, 
and therefore obtained the same scaling as we did in our quiet Sun derivation, which made the 
same assumption. 

The weaker dependence we find on photospheric flux concentration, Fh R 3 / 2 , may be in 
better agreement with observations. Some studies have synthesized EUV and X-ray images from 
coronal equilibrium fields with ad hoc heating fluxes Fh ~ B u [ 291 30]. These found the most 
reasonable matches with observations when v was near unity. 

Heat flux, Eq. (USD, also scales with the large-scale twist, a, in the AR, since that is assumed 
to be the source of current. This is a a natural dependence for any model since more twisted fields 
contain more magnetic energy to be released by reconnection. A statistical study of active regions 
by Fisher et al. m found no evidence for a scaling with the mean a of an AR. This could mean 
that reconnection is not the primary source of heat for ARs, or that the currents across which 
reconnection occurs are driven on scales smaller than the entire AR. 

Whether it supplies the majority of coronal heat or not, magnetic reconnection is clearly gen¬ 
erating heat in the corona at some rate. We have quantified this in Eq. © as the rate of electro¬ 
magnetic work attributable to topological change. This is not necessarily Joule heating. In several 
more detailed reconnection models this work is transferred primarily into bulk kinetic energy either 
in supersonic reconnection outflows, which generate heat through shocks [32] ]33l fl6 l or in MHD 
waves which ultimately damp to generate heat [33113‘5]. In view of the latter option, reconnection 
is not necessarily an alternative to heating by waves; rather it is a potential source of waves to heat 
the corona. 

This work was partially supported by a grant from the NSF/DOE plasma sciences program. We 
gratefully acknowledge two anonymous referees whose comments helped us improve the manuscript. 
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